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1 Introduction 

Observation of atmospheric |]l| and solar p[ neutrinos provide important indications that neu- 
trinos oscillate between different mass eigenstates. Interpreting each observation in terms of 
two-flavour mixing the oscillations of atmospheric neutrinos require Am^^^ ~ 2 x 10~^ eV^ and 
a nearly maximal mixing angle sin^ ^Oatm > 0.82. For solar neutrinos, several solutions of their 
deficit problem are possible. The vacuum oscillation (VO) solution requires Am'j.^i ~ 10~^° eV^ 
and sin^ 26^50^ > 0.67, whereas the the so-called MSW solution requires Am^^/ ~ O{10~^) eV^ 
and sin^ 29soi > 0.5 (large mixing angle solution - LAMSW) or 4 x 10~^ < sin^ 29soi < 1-3 x 10~^ 
(small mixing angle solution - SAMSW). Combining these results with non-observation of the 
dissappearance of z/g in the reactor experiments (in particular in CHOOZ [Q) sensitive to 
Am'^ > 10~^ eV^ and sin^ 26 > 0.2, the atmospheric and solar neutrino oscillations can be 
explained in terms of three known flavours of neutrinosQ the assumption which we adopt in 
this paper, as z/^ — >■ u^ and u^ — >■ u^ oscillations with u^ = ly^ if the mixing was two flavour one. 

With accumulating data, considerable attention has been focused on the determination of 
the full Maki-Nakagawa-Sakata (MNS) 3x3 unitary mixing matrix |^ which is conveniently 
parametrized by 3 angles (we set the phase 6 = 0) 

C1C3 S1C3 S3 

U = I -S1C2 - C1S2S3 C1C2 - S1S2S3 S2C3 I (1) 

S1S2 - C1C2S3 -C1S2 - S1C2S3 C2C3 

where Sj and q denote sin 6*4 and cos 6^4, respectively. 

Defining Am^ = m\—m{ with |Am^| = Aml^i and AM^ =m\ — m\ with |AM^| = Am^^^, 
the CHOOZ experiment measures (with AM^ ^ /S.rn? and sin^(Ar7i^L/4£') ^ 0) 

P{^e - ^e) = 1 - AUUl - Ul^ sin^ (^) (2) 

and the negative result of the search constraints S3 = f/13: S3 < 0.2. 

For the solar neutrino oscillations we have (averaging sin^(AM^L/4ii^) to 1/2): 

Pzy.zil'e -^ I'e) = clP2x2{l^e ^ I'e) + si (3) 

(for this formula to account also for the MSW solution one has to mutiply the electron density 
in the Sun in the 2-fiavour probability P2X2 by c\ |0|). In the limit S3 = 0, 6soi = Oi- However, 
since in general U31 7^ 0, P(z/e — > z/,-) 7^ 0. 



For the atmospheric neutrinos 



fAM'^L\ ^ ^ ^ ^ ( AM'^V 



P{v, - v^) = A{U,,U,,f sin^ y-^^j = Hclcl sin^ (^=^ ) (4) 



^We leave out the not yet confirmed LSND result jQ. All three measurements can be explained only if there 
exists sterile neutrinos, e.g. light right-handed ones 



so, in the limit S3 = 0, 9atm = ^2- 

Thus, up to some uncertainties in \Uis\ = {s^l, for each mentioned earher solution to the 
atmospheric and solar neutrino problems we can infer the gross pattern of the 3x3 mixing 
matrix. The striking features are large mixing angles 62 and, for VO and LAMSW solutions, 
61. Measuring S3 is very important, as the dependence of the fitted Si and S2 on S3 is not totally 
negligible [g g |Tg. 

Large number of papers addressed the issue of theoretical explanation of neutrino masses 



and mixings and of incorporating it into some global solution to the flavour problem |Tl|, [1^, 



|T^, |T^, |T5|, |T6|, |T^ |T^, ^, [2^, ^, ^ . One attractive possibility for giving the three flavours of 
neutrinos small Majorana masses is the seesaw mechanism which generates in the effective low 
energy theory (MSSM or the SM) a dimension five operator. In the two-component notation 
it reads: 

Cd^m 5 = -^C^\Hn{m') + h.C. (5) 

where a is generation index, I is the ordinary lepton SU{2) doublet, H is the hypercharge 
+1/2 Higgs doublet (i.e. if'-^-* in the MSSM) and the matrix C is dimensionless because we 
have factorized out the overall scale M of the heavy, right-handed neutrino mass matrix. This 
operator, after the electroweak symmetry breaking, is the origin of the left-handed neutrino 
Majorana mass term 

jC, ^ass = ~m'i'i^''i^' + h.C. (6) 

with m^ = [v'^ /AM)C where v/ \/2 = {H2). After the unitary rotation z/" -^ V^^z/^ diagonalizing 
the neutrino^ mass matrix C 

V[CVl = Cd = diag(Ci, C2, Cs) ~ diag(mi, ma, ms) (7) 

the charged current weak neutrino interactions depend on the MNS unitary matrix (||), U = 
E\Vl and the rotations El are defined by 

eImIEl = Hi ^ diag(/^,^ h% hi) (8) 

{Ail = Y^Ye is the square of the lepton Yukawa coupling matrix and hi are its eigenvalues). 
It is this matrix whose elements are probed in the neutrino oscilation experiments. In the 
effective theory below the scale M, with the right-handed neutrinos decoupled we can work in 
the charged lepton mass eigenstate basis, fixed by the diagonalization of the charged lepton 
Yukawa matrix at the scale M and then U = Vl. 

The measured neutrino mass matrix 171,^ = [v"^ / M)U C dU'^ is linked to the fundamental 
mass generation mechanism by two steps. The first one (in the top-down approach) is the 
renormalization group (RG) evolution of the neutrino and charged lepton Yukawa coupling 

^Wc will always assume that 7713 is positive and allow for negative signs of rai and/or rn2. 



matrices Y^ and Y^,, respectively, from the scale at which they are fixed by some theory of flavour 
(presumably at the GUT scale) down to the right-handed Majorana mass scale M < Mqut- 
At that scale the right-handed neutrinos with the Majorana mass matrix M are integrated out 
and one obtains the Standard Model (SM) or its supersymmetric extension (MSSM) with the 
Majorana mass matrix of the light neutrinos given by the operator (J^). The effective neutrino 
Majorana mass matrix C{M) at the scale M is dependent on both, neutrino and charged lepton 
original Yukawa textures at Mqut and on M. 

The second step consists of the renormalization group evolution of this operator down to 
the electroweak scale Mz , which provides unambiguous mapping of the pattern at the scale M 
into the measured pattern at Mz scale. Once we choose to work with the SM or the MSSM 
the physics below the scale M is almost unambiguous^] and it is interesting to study the effects 
of quantum corrections to the neutrino mass matrix C{M) summarized in its renormalization 
group (RG) evolution down to Mz- We find it convenient to work directly with rria's and f/a^'s 
rather than with elements of the matrix m^j since, as we shall see, the evolution of mass 
eigenstates and mixing angles often allows for easy qualitative discussion. 



In the present paper we extend earlier discussions |^, [2^, ^ by, first, writting down the 
renormalization group equations for the neutrino mass eigenvalues and mixing angles between 
the three fiavours of neutrinos, in the MSSM (and the SM). We point out the main differences 
between the evolution with 3x3 mixing or 2 x 2 mixing. Next, we discuss in some detail 
the solutions to the RGE's for mixing angles and their dependence on the mass eigenvalue 
patterns: m\ ~ AM^ 3> m\., mf (hierarchical), ml ^ m\ ~ AM^ ^ m\ (inversely hierarchical) 
and mf ~ m\ ~ m\ ~ (9(AM^) or larger (degenerate), which are consistent with the measured 
mass squared differences. We identify several infrared (IR) fixed points in the RG equations 
for mixing angles, which can be reached for some range of neutrino masses and are compatible 
with the present experimental information on the mixing angles. 

Finally, we discuss the stability with respect to quantum corrections of the three neutrino 
mass eigenvalue patterns which turns out to be correlated with the character of the evolution 
of the angles. In our formalism, we can reproduce in a very simple way the earlier results 
[p^ , p5| , p6| , in particular, those on the possibility of the VO solution with the stability of the 
inversely hierarchical or degenerate patterns. 



2 RGE for neutrino masses and mixing angles 

The importance of the scale dependence of the coefficient C of the operator (1^) have been first 



emphasized in ref. [27, E8|. The derived there renormalization group equation in the MSSM 



^In the of the SM there is a shght dependence of the overall scale of the neutrino masses on the Higgs boson 
quartic coupling A (i.e. on the Higgs boson mass) and in the MSSM on the assumed scale of supersymmetry 
breaking which we will ignore in the discussion which follows. 

"^For similar approach in the quark sector see P3|. 



real 



dsQ 



jC = -KC - CMl - {MlfC (9) 



with A^g = H^ (diagonal) in the basis we are working and where 
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K 



-^gf-Ggl + QTiMl 



(10) 



Aif^ — '^u'^u is the square of the up-type quarks Yukawa couphng matrix, gf = {5/3)gy, 
t = (I/IGtt^) log{M/Q) and M is the large Majorana scale. In the SM the correct form of the 
equation was given in ref. PU| and has the same form with —{6/5)gl -^ +2A (A is the scalar 



quartic coupling), —6g2 -^ — Sfifl, GTiM^ — > Tt{6M^ + 6M'^ + 2Ml) and the coefficients of 
the last two terms are +1/2 (instead of —1). 

Equation (H) can be elegantly solved ^T], since in the absence of right-handed neutrinos 
below the scale M the matrix El defined in (|]) does not run. Thus, in the basis in which the 
leptonic Yukawa matrix is diagonal the equation (P) has the obvious solution 

C{t) = IkIC{0)I (11) 

where X = diag(/e, /^, It), and 

Ik = exp (^- 1^ K{t')dt'^ , h = exp (^- J^ hl{t')dt'^ . (12) 

The only role of the factor Ik is to change the overall scale of the neutrino masses during the 
evolution. For M = 10^'^^^^ GeV, Ik ~ 0.9 — 0.6 with smaller values for lower tan/? due to the 
enhancement of the top quark Yukawa coupling in K. 

Although the solution to the RG equation for the matrix C is simple, qualitative features 
of the running of the mass egenvalues rria and the MNS mixing matrix U are often masked by 
the diagonalization procedure. In this paper we derive the RGE directly for rua and U. 

Following the method of ref. [^ , the RG equation for the matrix Vl = U (in the charged 
lepton mass eigenstate basis) can be written as 

—U = —Ue„ with eI = —e^ (13) 

where the matrix e^ is antihermitean, to preserve unitarity of U . It is determined by the 
requirement that the matrix C{t) is diagonalized by U{t) at any scale and the RHS of 

j^iU^CU) = jCd = -sICd - Cne. - KCd - CdU^HIU - U^HIWCd (14) 



^The MSSM RGEs given in ref. |g^ allow to treat also the case in which squarks and/or gluino are much 
heavier than sleptons, charginos and neutralinos so that the decoupling procedure ||29(| can be employed; in this 
case there are four different (in component fields) operators which mix below the squark/gluino treshold. Above 
it one has (in the notation of ref. [|8| ) cf' = 2c°2 = 2c2i = cf' = C""^ and the four equations of ref. ||2^ merge 
into the one quoted here (derived independently also in ref. pQ|). 



is diagonal. For real matrices C and Ai^ which we consider here, the MNS matrix U is real and 
orthogonal and e^ is antisymmetric. We get therefore 



.ab ^0. + (-^h ijjT tt2i 



C = -7^-^iU' H'eU)at for a ^6, 5^ = 0. (15) 

^a — ^b 



The running of the eigenvalues is then given by 



^Ca = - E (^ + '^hlU^a)Ca, a =1,2,3 (16) 

b=e,fiT 

and the running of elements of the MNS matrix is given by 

jU.^ = T.^^UUU^H!UU. (17) 

at ^^^ Cd - Cb 

If two of the three eigenvalues, say Ca and Cb, are equal at some scale t there is the obvious 
freedom in choosing the matrix U{t), corresponding to the redefinition U{t) -^ U{t) = U{t)R 
where i? is a rotation in the ab plane. This freedom is eliminated by quantum corrections, since 
for the evolution, R has to be fixed by the condition 

{U^HlU)ab{t) = {U^HlU)ba{t) = (18) 

so that equation (|1^) in nonsingular. 

Eqs. (|1^|1^) give directly the running of the the measurable parameters. It is straighforward 
to derive the equations for the three independent parameters Si, $2 and S3. Neglecting hf, and 
h^ Yukawa couplings we get: 

Si = -Ci(siS2 -CiC2S3)(-CiS2 -SiC2S3)y42l/l^ (19) 

- SiCiC2S3(siS2 - CiC2S3)v43i/l^ + clc2S'i{-CiS2 - 510253)^32/1^ 

S2 = Sicl{siS2- CiC2S^)A-iihl- Cicl{-CiS2- SiC2S^)A^2hl (20) 

S3 = -CiC2C3(SiS2 -CiC2S3)A3i/l^-SiC2C3(-CiS2 -5102^3)^32/1^ (21) 

where Aab = {Ca + Cb)/{Ca - Cb) = [rria + mb)/{ma - rrib). 

Eqs. ( P!B| - PT| ) give us several immediate results. Neglecting the electron and muon Yukawa 
couplings, the solutions for squared mass eigenvalues read: 

ml{t) = mlWl exp (- |J Ahl{t')UUt')dt^ (22) 

Observe that, since I\ <1, the masses always decrease top-down. We also see that the possi- 
bility of some change caused by the evolution in mass pattern resides solely in the differences 
in the mixing matrix elements U^a and their RG running. With h"^ ^ (tan/5/100)^, t ^ 0.12 for 



M = 10^° GeV and U^^ typically varying between and 1/4 (except for U^^) the exponent is 
at most of order oi e = hi log{M/Mz)/l6n'^ ^ tan^^g x 10"^ < 2.5 x lO'^ for tan/3 < 50. We 
can then estimate the changes in the mass squared differences: 

Amlit) ^ mlit) - ml{t) = Aml{0) - {Vaml{0) - Vtmlme (23) 

where we have neglected Ik which is always close to 1 and the factors rja^b) > are typically 
in the range — 2, depending on the values of U^a factors and their evolution. Taking (for 
definiteness) Am^j,(0) = 0, we see that the evolution of Am'^f,{t) is limited by m^(0)e or mf (0)e, 
i.e. by the value of the larger mass. 

A spectacular feature of eq. (|l3) or eqs. (p!9| - pT1) is the existence of fixed points at Uab = 0. 
The character of the fixed points (UV or IR) and the rate of approaching them depend on the 
mass factors Aab- This point will be discussed in detail later on. We note the possibility of a 
large factor Aab, for nearly degenerate and of the same sign eigenvalues rria and rub, which can 
permit a rapid approach to fixed points and drastically change the pattern of mixing at Mz 
during the evolution from M. f\ In the solutions ([11]) for the matrix C this point is encoded in 
the diagonalization procedure one has to perform after the running. The equivalence of the two 
approaches can be proven by diagonalizing the matrix C{t + At) using the ordinary perturbation 
calculus. In the first order in At one obtains then precisely eq. (p!?!). This approach also can 



serve to justify eq. (jJaj} . 



Finally we comment on two flavour mixing (which is sometimes good approximation to the 
more complex three flavour mixing). The 2x2 MNS matrix is an ordinary orthogonal rotation 
matrix. The general formalism applied to the mixing between say, second and third generation, 
gives 

j- sin^ = sint^cos^^ ^' ^ ^' {hi - hj). (24) 

at C3 — C2 

It is then straighforward to write down the equation for the evolution of sin^ 2'd: 

4-sin2 2^ = 2 sin' 27? cos 2^ ^'^^' (/i^ - hj). (25) 



Observe also that eq. (|25D, contrary to the original one (0), cannot be used for the initial 
condition cos2'(9 = 0. This is because cos2'(? = is the point at which the uniqueness of the 
solution of the differential equation ( p5D is violated. It has there two solutions: a t-dependent 
one which is the solution also to eq. (|2^ ) and the second one, sin' 2^ = 1, which does not 
solve eq. (U^j. Thus, the previous claims that the maximal mixing sin' 26' = 1 is stable against 



the RG running are not correct [pQ, 33, 34, 21, 35]. It is easy to see that, for the two-flavour 



evolution to be good approximation to the full three-flavour evolution of 62 as given by eq. 
(pOD, two conditions must be met: i) 63 has to be small and ii) the factors A^i and A32 have to 
be comparable. 



^The possibility of strong changes in the mixing pattern due to the evolution has been pointed out earlier 
for 2x2 mixing $^. 



3 Evolution of mixing angles 

The equations derived in the previous section can be used for a systematic analysis of the RG 
evolution of the mass and mixing patterns which are consistent with neutrino experimental 
data. To fix the framework we universally impose the following constraints at the electroweak 
scale: AM^ = 10"^ eV^, sin^ 26^2 > 0.82 and sin^ 9^ = sj < 0.2. Furthemore, for each solution - 
VO, LAMSW and SAMSW - to the solar neutrino deficit problem we assume the corresponding 
value of sin^ 29i and consider all mass configurations that are consistent with the respective 
Am^. The list of acceptable low energy mass patterns is, therefore, as follows: 

I AM^ ^ m\ 3> m\,m\ (hierarchical) with two distinct subpatterns: 

Am? K. m\,^^ ^ "^1(2) °^ 

7722(1) ~ ^1(2) with TTigM) ranging from m^M^ > Am^ to 7712(1) ^ Atti^ (the latter being, 
strictly speaking, possible only in the VO case; for SAMSW and LAMSW solutions, 
with our representative numbers A7rt^ = 5 x 10~^ and 5 x 10~^ eV^, respectively, 
777,2(1) -^ Am? violates the inequality AM^ ^ 77i| ^ ttIi, m\ and one gets the pattern 
III below). It is also necessary to distinguish the cases mim,2 < and mim2 > 



II AM^ ^ m\^m\ ^ 77I3 (inversely hierarchical). Here, also, the evolution depends on the 
choice 777i7772 < (witli Atti^ > or < 0) or mim2 > (with Am'^ > or < 0) 

III TTig ^ m\ K, ml and all of the order or larger than AM^ (degenerate). Here the evo- 
lution depends on the relative magnitudes (four different orderings satisfy our general 
constraints) and relative signs of the masses (we choose 7713 > and for each ordering 
there are four inequivalent sign combinations). 

Since the evolution of angles is less dependent on the evolution of mass eigenvalues than 
the other way around, it is convenient to begin our discussion with the former. It can be 
classified into several universal types of behaviour, depending on the magnitude of the factors 
Aab in eqs. (|T9|-^T]). We note that, neglecting the effects of mass evolution, all possible mass 



configurations contained in patterns I— HI give one of the following four structures: 

a) A31 ^ A32 and IA31I ^ |v42i| ^ 1 

b) A31 ^ A32 and IA21I > IA31I, \A2i\ > 1 

C) ^32^ ^21^0, 1^31 1 >1 
d) A3i^A2i^0, 1^32 1 >1 



For a) and 6), since A^i ^ A^2, equations ([T9| - pID reduce to: 

51 = -Cl{SiS2-CiC2Ss){-CiS2- SiC2S3)A2lhl-CiS2C2S3A32hl, 

52 = S2clA^2hl, (26) 

2 2 /I 7,2 



Thus, the evolution of the atmospheric mixing angle is as in the two-generation equations 
or (^51). The first two patterns, hierarchical and inversely hierarchical give the Aab factors as 
in a) or h). In addition, since for patterns I and II IA32I ~ 1, the evolution of both, S2 and S3 
is very weak (note that S2 = and/or S3 = are fixed points). Denoting 

e. = exp {- J^ 2A;2{t')hlit')dt'^ ^ exp (-2^32(0)6) (27) 

the solution for sl(t) reads 

slit) = sl{0)/\sl{0) +4(0%] (28) 



and yields 



sm'2e2{t) = ^r sin' 2^2(0)/ s^iO) + 4(0% 



2 



^ sin' 2^2(0) [1 + 2 cos 2^2(0)^32(0)6] + 0{Al^e'^). (29) 

The effect of the running for sin' 29atm is a 2.5% change for extreme value of tan/5 ~ 50. Since 
there is no large factor involved, the logarithmic approximation (the last line in eq. (|29D ) is 
very good. This is important for the evolution of the masses, discussed in the next section. 

Moreover, we observe that, in the very good approximation of constant C2, the solution for 
S3 is also of the form (p8|), with S2(c2) -^ ^3(03) and ^32 — > 02^32 in the definition of ^^. 

For the solar mixing angle, in the approximation S3 = 0, we have 

Sl = siclslA2ihl (30) 

and the solution for si is of the form ( P5| ) too, with S2(c2) -^ si(ci) and 

e. - e = exp (- 1^ 2sl{t')A2^{t')hl{t')dt'^ . (31) 

With further specification to the structure a), occuring for Am' ~ "^2(1) -^ "^1(2) ^^'^' i^ 
general, for mim2 < 0, the evolution of the angle Si is very weak, independently of the chosen 
solar solution. Since all \Aab\ ^ 1 the dependence on S3 is weak (for s| < 0.2). We conclude 
that in case a) the evolution of all mixing angles is negligible and the existence of the fixed 
points at Uab = in eqs. ([T9| - |2l|) is irrelevant. 



The evolution of si is dramatically different for \A2ie\ ^ 1, i.e. for m^^-. > |Am'| and 
17111712 > 0. The presence of the fixed points becomes relevant but, since the evolution of S2 
and S3 is still very weak (IA32I ~ 1), they manifest themselves as two approximate fixed points 
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A2i£(M) 



A2i£(M) 



Figure 1: Approach to the IR fixed point for the angle 6i as a function of A2ie{M) for two 
different values of S3: a) s| = 0.2 and the hierarchy Imsl ^ |m2| Z \mi\, b) S3 = (0.025)^ and 
the hierarchy \m2\ ~ \Tni\ ^ I'tt-sI- In both cases mim2 > and we have taken sin^ 2^2 = 0.9 so 
that in case a) fixed point corresponds to sin^ 29i = 0.8 and in case b) to sin^ 29i = 5 x 10^'^. 
Different lines correspond to different choices of Oi{M) As explained in the text, for |v42i|e ^ 1 
and 1^432! ~ I ^31 1 ~ 1, the evolution of 62 and 63 is negligible. 



of the RG equations for Si. One can easily check (for instance, by considering the equation 
for dtanOi/dt = (l/c?)si) that for A21 > (i.e. for Am^ > 0) the point f/31 = is the UV 
fixed point and U32 = is the IR fixed point. For A21 < (i.e. for Am^ < 0) the situation is 
reversed. It is also interesting to notice that in the limit S3 = we can follow analytically the 
approach to the fixed points. In this limit the RG equation for si has again the form ( pOD with 
the solution of the form (P^ with the replacement (0). Therefore, for A21 > 0, in the top-down 
running the factor ^' — > exponentially with growing A2ih'^\og{M/Mz) and, consequently, we 
obtain Si(t) = ±1 (depending on its initial sign) and approach IR fixed point at t/32 = 0. In the 
bottom-up evolution we approach sf(t) ~ exponentially, i.e. the UV fixed point at f/31 = 0. 
For A21 < we get the reversed situation, in accord with our general expectations. 



Several comments are in order here. First, it is useful to remember the relation between 



sin^ 29i and sin^ 2^2 as the function of S3 at U31 = and U32 = (i.e. for tan 9i 
and tan^i = — tan6'2/s3, respectively). For both we get the following: 



sin^ 2^1 



si sin^ 2^2 

{sWs + sly 



with 



^(l±Vl-sin^2^,i„ 



S3/ tan 6*2 



(32) 



Thus, both fixed points at IR are at present experimentally acceptable with the mixing close 

0.2 and with the SAMSW solution for S3 ~ 0. This is interesting in the 



to bi-maximal for S3 ~ 
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sense of the independence of initial conditions at the scale M. 

Secondly, it is interesting to estimate the values of A21 and tan (3, for which the approach 
to the IR fixed points is seen. This is shown in Fig. |I|. We can estimate that for approaching 
the fixed point during the evolution in the range (M^, M) one needs A2ie{M) > 3, i.e. for 
tan/? = 20 one needs mi p^ m2 ^ 10~^ eV for the VO and rrii ^ 777-2 ~ 0.01 eV for LAMSW or 
SAMSW solution. 

Finally we note, that from the point of view of the initial conditions at the scale M, the UV 
fixed point looks not realistic as the neglected muon Yukawa coupling /7^ quickly destabilizes 
it during the evolution. We conclude that for hierarchical and inversely hierarchical mass 
patterns the evolution of the mixing angles is either very mild or shows (for |A2i|e ^ 1) fixed 
point behaviour. With higher precision experiments, it will be very interesting to confirm or 
disprove the IR fixed point relation between the angles. 

The evolution of mixing angles in the degenerate case, m^ ~ tti^ ~ ml ~ (9(AM^) or larger, 
partly falls into the same classes of behaviour. Indeed, as long as ^31 ^ A32 and 1^21 1 ^ I ^31 1 5 
I ^21 1 3> 1, the angles evolve according to the same equations (p6|). One can easily identify the 
eight mass patterns of the degenerate case that fall into this cathegory: the sufficient condition 
is that mi and 7772 are of the same sign. For the evolutions of si we then closely follow the 
two possibilities (depending on the sign of A21) discussed for the first two hierarchies with 
777i7772 > 0. Wc siuiply uotc that larger values of |^2i| are generic for the present case and, 
as seen in Fig. |I|, the approach to the fixed points is faster. However, the evolution of S2 and 
S3 are mild only if mi and 777.2 are negative (^431 ^ ^132 ~ 0). For positive mi and 7^2, we 
have A^ie ^ ^4326 and |v432e| can be much larger than 1 (depending on the overall mass scale 
0.1 eV < mi < 2 eV) so that the evolution of S2 and S3 is no longer negligible. According to 
the solution (pS]), S2 is exponentially focused to S2(t) = or S2(t) = ±1, depending on the sign 
of ^32, and on the direction of the evolution but independently of the values of Si and S3. We 
stress again that S2(0) = ±l/-\/2 i.e. siii^ 29atm = 1 is not stable. The angle S3 behaves in a 
similar way. Both angles can approach the value for which sin^ 29i = (7 = 2, 3), corresponding 
to the fixed points of their respective equations. This is experimentally acceptable as IR and 
UV fixed point for S3 but only as UV for S2. Since UV fixed points are unstable and small 
corrections push the evolution towards the IR ones, we conclude that the pattern ^31 ^ A32, 
\A3i\ ^ 1 is unacceptable in the regime in which the approach to the fixed points is relevant. 

The remaining degenerate mass patterns can be classified according to the relations A32 ~ 
^21 ~ or ^31 ?» A21 ~ 0. Consider first A21 ~ A32 ~ 0. The equations governing the 
evolution of the mixing angles can be approximated as 

Si = -SiCiC2S3(siS2 - CiC2S3)A3i/l^, 

52 = SiC2(SiS2 -CiC2S3)A3i/l^, (33) 

53 = -ClC2C3(siS2 - CiC2S3)y43i/7^. 

These equations exhibit IR quasi-fixed point behaviour for yl3i < corresponding to U31 = 0. 
As before, at the fixed point the angles satisfy the relation S3 = tan^itan6'2. We recall that 
this relation is consistent with present experimental information. Since Si is proportional to Si 
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-A3i£(M) 



b) 



sIn02(M) 



Figure 2: a) Approach to the IR fixed point for the degenerate mass pattern \m2\ ~ \mi\ > 
Imsl and nii > 0, 1712 < 0, as a function of — y43ie(M) for initial values of S2(M), s^i^M) 
and the masses such that Am'^{Mz) = 5 x 10"^ eV^, AM'^{Mz) = lO'^ eV^ s|(Mz) = 
0.2 and sin2 2^i(Mz) = 0.85 (LAMSW); at the fixed point sm'^2e2iMz) = 0.9. Different 
lines correspond to different choices of 62{M) b) Correlation of S2{M) and s^lM) which give 
experimentally acceptable angles at the fixed point. Parameters at Mz are as in panel a). 

and suppressed by S3, the running of si is weak. The IR fixed point is reached due to strong 
running of S2 and S3. The rate of approaching the IR fixed point by S2 is illustrated in Fig. |^a. 
For Asie 3> 1, S2{Mz) is strongly focused at ±1. Thus, mass and tan/5 configurations leading 
to A^ie ^ 1 are unacceptable. For ^432 < we get IR fixed point in U32 = 0. 

In summary, with all \Aab\ ~ 1 the evolution of the mixings is negligible. For \A2i\ ^ 1 
and 1^31 1, 1^32! < 1, or for ^32(^31) <^ — 1 and ^431(^432) ~ A21 ~ the infrared fixed points 
are reached during the evolution, independently of further details of the mass matrices. We 
also note that for |y42i| ^ 1 only si runs to assure the fixed point relation, so already the 
initial values for S2 and S3 at the scale M have to be close to their experimental values. For 
^31 ^ —1 or y432 <^ —1, S2 and S3 evolve strongly and the evolution of Si is weak. To assure 
consistency with experimental data, the initial value of Si has to be close to the experimental 
value and the initial values of the other two angles have to satisfy certain relation (see Fig. ^). 
Special examples of the IR fixed point behaviour are the textures with exact degeneracy of 
some masses, which assure the IR fixed point relations already at the scale M. They will be 
recalled in the last part of the paper. 



11 



4 Evolution of mass eigenvalues 

The guiding principle for understanding the effects of evolution of mass eigenvalues is eq. (pSf ) 
and the evolution of the angles. It follows from eq. ( P^D that quantum corrections | Am^^(M^) — 
Am^j(M)| are limited from above by 0{m'^,^^e) where e ~ tan^/3 x 10~^. In fact, generically 
this upper bound is saturated (in the sense of order of magnitude) for those mass patterns 
which give fixed point solutions for mixing angles (at all scales or reachable after the evolution) 
since then in general (f/|i) 7^ (^^32) (where (. . .) means average over the evolution). Quantum 
corrections orders of magnitude smaller than this upper bound can be naturally guaranteed 
only for weakly evolving angles and with [/|^ cs U^^. Those two rules determine the order of 
magnitude of quantum corrections to any mass pattern of interest. We illustrate them with 
several examples. 

The fully hierarchical pattern, with AM^ ^ m| ^ |A??i^| ^ ^^1(2) ^ "^2(1)' ^^ obviously 
not altered by quantum corrections. Similar conclusion holds for ml ^ m\ ~ (9(A??i^). More 
interesting effects may appear for the hierarchy AM^ ^ m\ ^ in\,m\ ^ Am^ which is 
possible for the VO solution. It is then interesting to ask the same question which is usually 



asked ^, 0, ^, 0, ^ for the inversely hierarchical and degenerate patterns: is the initial 
condition irti = ±m2 at the scale M compatible with the measured (at the scale Mz) Am^ ? 
We consider first mi = m2 which places us in the structure b) considered previously for the 
evolution of the angles. We note that eq. ([TsD imposes (for h^ = he = 0) at the scale M the 



fixed point values U^i = or f/32 = 0. It follows then from the evolution of the masses, eq. (^2]) 
and our discussion of the fixed point solutions, that both choices are IR fixed points, so remain 
stable during the evolution. For f/3i(32) = we get Am|]^(M^) = ±mw2-)(M)4e?732/3]^), where 
^32(31) ~ ^(-^i)- SO' fo^ "^1 — "^2 ~ lO^*^^^^) eV and tan/3 < 2.5 one obtains the right order of 



magnitude for the VO solution |25 . 



Choosing nii = —m,2, the angles remain uncorrelated with the masses (structure a)) and 
stable. For m^f.^) ~ 10~^^~^^ eV^ it is easy to choose angles so that Am|(M^) ~ (^(emwg)) ~ 
10^^*^ eV^ in a large range of tan/5 values. 

One should also note that for the hierarchy I with m,l{M) = m2{M) getting the right Am'^ 
for SAMSW or LAMSW solutions is only marginally possible for m^ > 10^^ eV^ and large 
values of tan/3 (> 50). 

We consider now the inverse hierarchy AM^ ~ m^ ^ m\ ^ m\. For mi{M) = m,2{M) we 
can repeat the previous discussion remembering, however, that now m^/g) ~ 10^'^ eV^. Thus, 
with such initial conditions we cannot reproduce Am"^ of the VO solution but we naturally get 
Am^ of the SAMSW (for 2.5 < tan/3 < 25) or LAMSW (for 5 < tan/3) solutions [||]. And, of 
course, the mixing angles satisfy one of the fixed point relations for all scales. Note, however, 
that as stressed earlier the exact degeneracy at M is not necessary to reach the IR fixed point 
at Mz 

The choice mi = —m,2 guarantees mild evolution of the elements U^a, Uzaif) = ^3a(0)+C(e). 
Since m\ ~ m,\ ~ 10^^ eV^ and Am\^{Mz) ^ 0{em,\,2){^)) ^^ note that choosing the angles 
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in agreement with experimental information we can reproduce the right order of magnitude 
for Am^ for the LAMSW or SAMSW solutions (for roughly the same ranges of tan/5 as for 
mi = 1712). On the other hand, for the VO solution, one needs Am^iiMz) ~ (9(e^m^/2)(^)) 
[p^j . The absence of C(emw2)(M)) contribution to Am? requires U^i = U^2 which has two 



solutions: 



.3 = -(^)^^-^tank-^) or (34) 

\C2/ Ci + Si C2 V 4/ 



,3=|£a)£l±il._£icot(9,-f) (35) 

\C2/ Ci - Si C2 V 4/ 



independently of tan/?. 



The easiest way to find consistently the mass splitting in order e^ is to solve the coupled 
equations (pJ| , p^pl| ) in order e^ by using the so-called Banach's principle i.e. by substituting 



in the RHSs of the RGEs their solutions obtained in the order e and integrating once more. 



One then gets (for definiteness we assume (0) and following ref. p4| take the neutrino mass 
spectrum in the form ~ (—1, 1, 2;)): 



m^ = ml{l-AUle + 8U^ie'-AUlclcl^e' 

,2/1 ,irr2 , I orr4 ,2 ^7-72 2 2^ + ^,2 



m'2 = mtil-AUiie + 8Ulie'-4Uii44—-e'). (36) 



z 



where U^i = U^2 = (■S2/(si + Ci))^. We have used the approximation of constant h^. and 
absorbed the overall renormalization introduced by the factor K in eq. (|16D in m^. We obtain, 
therefore 



,2 _ ig_,2„2„2^ '^3 \ ^ ^2 

Xl + SiJ Z^ — \ 



Am' = -lQmtsi4 ^- — — -e^ (37) 



For S3 = and 31 = 01 = 1/2 this coincides with the result derived in ref. p4|. Since 
AM2 = m'^{z^ - 1) + C(e) we get H 



^-'--''[Ai^)-^'^^'^[wt7i)- ^''^ 

We note that for inversely hierarchical pattern, z = 0, 0{e^) contribution to Am"^ vanishes 
too. For z ~ 0{1) the above considerations are applicable also to the degenerate pattern III 
provided |v43i(2)|e < 1 (see Fig. ||a). With e < 2.5 x 10"^ for tan/5 < 50 we get ^31(2) < 50 and 
z < 49/51. For a large range of parameters it is easy to obtain Am^ consistent with the VO 
solution. For the degenerate pattern III with ^31(2)6 ^ — 1 the expansion in e (in fact for the 
mixing angles this is an expansion in £^432) breaks down and this mechanism cannot work. 

In general, with the degenerate pattern III, the potentially relevant simple textures are 
m^ > \m,i\ = \m,2\ or |?7ii| = \m,2\ > m,^. 
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The exact degeneracy of absolute values {m^l = \m2\ = \mi\ would not result in two distinctly 
different AM^ and Am^ after evolution. With mi and m2 both negative we are in the IR fixed 
points and encounter similar situation as for the inversely hierarchical case with mi = m2- The 
result Amii(Mz) ~ 0{tm\^^-^{M)) is compatible with the LAMSW and SAMSW solutions but 
not with the VO solution. However, the relation Am\i{Mz) ~ Amlgl^z) ~ AM^ can be 
simultaneously satisfied only if all masses squared are of order AM^ ~ lO'*^^"^-* and therefore 
we need 2.5 < tan/3 < 15 (5 < tan/3) for the SAMSW (LAMSW) solution. As explained in the 
previous section, both masses positive are unacceptable in the regime in which the fixed points 
are relevant. 

For nil = —m2 we get ^421 = 0, A32 ~ for positive nii and A21 = 0, ^31 ^ for negative 
nil, with 1^31 1 ^ 1 and IA32I ^ 1, respectively. Such initial conditions assure the reaching of 
IR fixed points at U31 = for Asie ^ — 1 or U32 = for A32e ^ — 1 and, from the point of view 
of quantum corrections to the masses, we are back to the same situation as for mi and m2 both 
negative. For |y43i(2)|e < 1 but ^31(2) <^ —1 the angles evolve very mildly and the previously 
described mechanism of stabilization can assure right Am'^ for the VO solution 



5 Conclusions 

In this paper we have derived RG equations directly for the mixing angles and mass eigenval- 
ues. These equations allow for easy qualitative discussion of the evolution of masses and mixing 



angles, which systematizes existing results [B^, 25, 24, 37, 3S, M and allows for their generaliza- 



tion. The equations for the mixing angles have IR fixed points at t/31 = or U32 = 0, which are 
reachable for several mass hierarchies consistent with the measured mass squared differences 
and characterized by one large factor \Aab\ = K^^ + mh)/{ma ~ m^)] ^ 1 with all other IAq^I's 
< 0{1). Both fixed points give relation (|3^) which is at present acceptable experimentally. Its 
further verification, by clarifying the LSND puzzle, measuring S3 and discriminating between 
large and small solar angle solutions is of obvious interest. 

Special examples of the IR fixed points behaviour are the textures with exact degeneracy of 
some masses. The mass squared differences consistent with the measured values can be obtained 
for the following textures: i) hierarchical with Im^] ~ 10~^'^ eV, mi = m2 and \mi\ ~ 10~*^^~^) 
eV (for VO), ii) inversely hierarchical with rrii = m2 ~ 10^^'^ eV ^ m^, (for LAMSW and 
SAMSW), in) degenerate with |mj| > (9(0.1) eV {i = 1,2,3) and nii = m2 < or nii = —m2 
with ^316 < -1 (for LAMSW and SAMSW). 

With the fully hierarchical masses {m^l ~ 10^^'^ eV, m"^ ~ Am"^ ^ nif and for the inverse 
hierarchy with mi = —m2 the evolution of the angles is very mild. The proper order of 
magnitude for Am"^ is possible for the LAMSW and SAMSW solutions. Similarly mild evolution 
occurs for the degenerate case with mi = —m2, 1^31(2) |e ^ 1, ^31(2) -C —1 which is consistent 
with the stabilization mechanism leading to Am"^ of the VO solution. 

We focused our discussion on the MSSM. The corresponding RG equation in the SM are 
obtainable by the replacements mentioned in the text. The UV and IR fixed points are in- 
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terchanged and the running of the angles is weaker by factor 2 than for tan/? = 1 in the 
MSSM. 
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